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Abstract. The nonperturbative renormalization group flow of Quan- 
tum Einstein Gravity (QEG) is reviewed. It is argued that there could 
be strong renormalization effects at large distances, in particular a scale 
dependent Newton constant, which mimic the presence of dark matter 
at galactic and cosmological scales. 



1 Introduction 

By now it appears increasingly likely that Quantum Einstein Gravity (QEG), the 
quantum field theory of gravity whose underlying degrees of freedom are those of 
the spacetime metric, can be defined nonperturbatively as a fundamental, "asymp- 
totically safe" theory (Lauscher 2002). By definition, its bare action is given by 
a non-Gaussian renormalization group (RG) fixed point. In the framework of the 
"effective average action" a suitable fixed point is known to exist within certain 
approximations. They suggest that the fixed point should also exist in the exact 
theory implying its nonperturbative renormalizability. 

The general picture regarding the RG behavior of QEG as it has emerged so far 
points towards a certain analogy between QEG and non-Abelian Yang-Mills the- 
ories, Quantum Chromo-Dynamics (QCD) say. For example, like the Yang-Mills 
coupling constant, the running Newton constant G = G{k) is an asymptotically 
free coupling, it vanishes in the ultraviolet (UV), i. e. when the typical momentum 
scale k becomes large. In QCD the realm of asymptotic freedom is realized for 
momenta k larger than the mass scale Aqcd which is induced dynamically. In 
QEG the analogous role is played by the Planck mass mp\. It delimits the asymp- 
totic scaling region towards the infrared (IR). For k ^> mpj the RG flow is well 
described by its linearization about the non-Gaussian fixed point. Both in QCD 
and QEG simple local approximations (truncations) of the running Wilsonian ac- 
tion (effective average action) are sufficient above Aqcd and mpi, respectively. 
However, as the scale k approaches Aqcd or mpi from above, many complicated, 
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typically nonlocal terms are generated in the effective action. In fact, in the IR, 
strong renormalization effects are to be expected because gauge (diffeomorphism) 
invariance leads to a massless excitation, the gluon (graviton), implying potential 
IR divergences which the RG flow must cure in a dynamical way. Because of the 
complexity of the corresponding flow equations it is extremely difficult to explore 
the RG flow of QCD or QEG in the IR, far below the UV scaling regime, by 
analytical methods. In QCD, lattice results and phenomenology suggest that the 
nonpcrturbative IR effects modify the classical Coulomb term by adding a confine- 
ment potential to it which increases (linearly) with distance: V(r) = —a/r + nr. 

The problem of the missing mass or "dark matter" is one of the most puzzling 
mysteries of modern astrophysics. It is an intriguing idea that the apparent mass 
discrepancy is not due to an unknown form of matter but rather indicates that 
we are using the wrong theory of gravity, Newton's law in the non-relativistic and 
General Relativity in the relativistic case. If one tries to explain the observed non- 
Keplerian rotation curves of galaxies or clusters in terms of a modified Newton 
law, a nonclassical term needs to be added to the 1 / r-potential whose relative im- 
portance grows with distance. In "MOND" , for instance, a point mass M produces 
the potential <fi(r) = —GM/r + \faa GM ln(r) and it is tempting to compare the 
ln(r)-term to the qualitative similar confinement potential in (quenched) QCD. It 
seems not unreasonable to speculate that the "confinement" potential in gravity 
is a quantum effect which results from the antiscreening character of quantum 
gravity (Lauscher 2002) in very much the same way as this happens in Yang-Mills 
theory. If so, the missing mass problem could get resolved in a very elegant man- 
ner without the need of introducing dark matter on an ad hoc basis. In (Reuter 
2004a, b) this idea has been explored within a semi-phenomenological analysis of 
the effective average action of quantum gravity. 

2 RG running of the gravitational parameters 

The effective average action Tk[g^ u ] is a "coarse grained" Wilsonian action func- 
tional which defines an effective field theory of gravity at the variable mass scale k. 
Roughly speaking, the solution to the effective Einstein equations STk/Sg^ = 
yields the metric averaged over a spacetime volume of linear extension k^ 1 . (From 
the technical point of view k is a IR cutoff introduced into the functional integral 
over the microscopic metric in such a way that only quantum fluctuations of wave- 
lengths smaller than are integrated out.) In a physical situation with a typical 
scale k, the effective field equation STk/Sg^ = "knows" about all quantum ef- 
fects relevant at this particular scale. For k fixed, the functional should be 
visualized as a point in the space of all action functionals. When the RG effects 
are "switched on", one obtains a curve in this space, the RG trajectory, which 
starts at the bare action S = Tk^oo and ends at the ordinary effective action 
r = Tfe^o- At the exact level, Tk contains all the infinitely many invariants one 
can construct from g^ u , their fc-dependent prefactors having the interpretation of 
scale dependent gravitational coupling constants. To become technically feasible 
most of the investigations using the effective average action formalism employ the 
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so-called Einstein-Hilbert approximation which retains only Newton's constant 
G(k) and the cosmological constant A(fc) as running parameters. If one intro- 
duces the dimensionless couplings g(k) = k 2 G(k) and A(fc) = A(k)/k 2 the RG 
equations governing their scale dependence read kdkg — @g{g, A), kdkX = f3\(g, A) 
with known beta-functions f3 g and The RG flow on the g-A-plane displays two 
fixed points: a Gaussian fixed point (GFP) at the origin, and the non-Gaussian 
fixed point (NGFP) at > 0, A* > which is necessary for asymptotic safety. 
The RG trajectories are classified as of Type la, Ha (separatrix), and Ilia de- 
pending on whether, when k is lowered, they run towards negative, vanishing, and 
positive values of the cosmological constant, respectively. In (Reuter 2004b) the 
very special trajectory which seems realized in Nature has been identified and its 
parameters were determined. This trajectory is of Type Ilia; see fig. ^ 
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Fig. 1. Nature's Type Ilia trajectory and the separatrix. The dashed line is a classical 
RG trajectory along which G(k), A(k) — const. (From (Reuter 2004b).) 

For k — > oo it starts infinitesimally close to the NGFP. Then, lowering k, the 
trajectory spirals about the NGFP and approaches the "separatrix", the distin- 
guished trajectory which ends at the GFP. It runs almost parallel to the separatrix 
for a very long "RG time"; only in the "very last moment" before reaching the 
GFP, at the turning point T, it gets driven away towards larger values of A. In 
fig. n the points Pi and P2 symbolize the beginning and the end of the regime 
in which classical general relativity is valid ("GR regime"). The classical regime 
starts soon after the turning point T which is passed at the scale fcx ~ 10 _30 mpi. 

In (Reuter 2004b) we argued that to the right of the point P2 there starts 
a regime of strong IR renormalization effects which might become visible at as- 
trophysical and cosmological length scales. In fact, within the Einstein-Hilbert 
approximation, trajectories of Type Ilia cannot be continued to the extreme IR 
[k — > 0). They terminate at a non-zero value of k as soon as the trajectory reaches 
A = 1/2. (Close to the question mark in fig. [I]) Before it starts becoming in- 
valid and has to be replaced by a more precise treatment, the Einstein-Hilbert 
approximation suggests that G will increase, while A decreases, as A /* 1/2. 
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The Type Ilia trajectory of QEG which Nature has selected is highly special 
in the following sense. It is fine-tuned in such a way that it gets extremely close to 
the GFP before "turning left" . The coordinates gr and At of the turning point are 
both very small: gr — At ~ 1CP 60 . The coupling g decreases from g(k) = 1CP 70 
at a typical terrestrial length scale of k^ 1 = 1 m to g(k) — 10~ 92 at the solar 
system scale of fc _1 = 1 AU, and finally reaches g(k) = 1CP 120 when k equals the 
present Hubble constant Hq. 

In fact, the Hubble parameter k — Hq is approximately the scale where the 
Einstcin-Hilbcrt trajectory becomes unreliable. The observations indicate that 
today the cosmological constant is of the order Hq. Interpreting this value as the 
running A(fc) at the scale k — Hq, the dimensionlcss A(fc), at this scale, is of the 
order unity: X(Hq) = K(Hq)/Hq = 0(1). So it is roughly near the present Hubble 
scale where the IR effects should have grown large. 

In principle it should be possible to work out the predictions of the theory 
for cosmological scales by an ab initio calculation within QEG. Unfortunately, 
because of the enormous technical complexity of the RG equations, this has not 
been possible in practice yet. In this situation one can adopt a phenomcnological 
strategy, however. One makes an ansatz for the RG trajectory which has the 
general features discussed above, derives its consequences, and confronts them 
with the observations. In this manner the observational data can be used in order 
to learn something about the RG trajectory in the nonperturbative regime which 
is inaccessible to an analytic treatment for the time being. Using this strategy, 
the cosmological consequences of a very simple scenario for the k — > behavior 
has been worked out; the assumption proposed in (Bonanno 2002) is that the 
IR effects lead to the formation of a second NGFP into which the RG trajectory 
gets attracted for k — > 0. This hypothesis leads to a phenomenologically viable 
late-time cosmology with a variety of rather attractive features. It predicts an 
accelerated expansion of the universe and explains, without any fine tuning, why 
the corresponding matter and vacuum energy densities are approximately equal. 

3 Galaxy rotation curves 

Given the encouraging results indicating that the IR effects are "at work" in cos- 
mology, by continuity, it seems plausible to suspect that somewhere between so- 
lar system and cosmological scales they should first become visible. In (Rcuter 
2004a, b) we therefore investigated the idea that they are responsible for the ob- 
served non-Keplerian galaxy rotation curves. The calculational scheme used there 
was a kind of "RG improvement" , the basic idea being that upon identifying the 
scale k with an appropriate geometric quantity comparatively simple (local) trun- 
cations effectively mimic much more complicated (nonlocal) terms in the effective 
action. Considering spherically symmetric, static model galaxies only, the scale k 
was taken to be the inverse of the radial proper distance which boils down to 1/r 
in leading order. Since the regime of galactic scales turned out to lie outside the 
domain of validity of the Einstein-Hilbcrt approximation the only practical option 
was to make an ansatz for the RG trajectory {G(k), A(fc), • • • } and to explore its 
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observable consequences. In particular a relationship between the fc-dependence 
of G and the rotation curve v(r) of the model galaxy has been derived. 

The idea was to start from the classical Einstein-Hilbert action and to pro- 
mote G and A to scalar fields:.? = Jd 4 x sj-g {R/G(x) - 2A(x)/G(x)}. Upon 
adding a matter contribution this action implies the modified Einstein equation 
G M „ = -A(x)g MV + SnG(x) (T^ + AT^) with AT MV = gL (D^D v -g^ D 2 ) G" 1 . 
In (Router 2004a) we analyzed the weak field, slow-motion approximation of this 
theory for a time-independent Newton constant G = G(x) and A = 0. In this 
(modified) Newtonian limit the equation of motion for massive test particles has 
the usual form, x(i) = — V</>, but the potential (f> obeys a modified Poisson equa- 
tion: 

V 2 </> = 47rGp cff where p cS = p + (8vr G)" 1 V 2 W. (3.1) 

Here it is assumed that describes pressureless dust of density p and that G(x) 
does not differ much from the constant G. Setting G(x) = G [l + AA(x)] we as- 
sumed that A/"(x) <C 1. Apart from the rest energy density p of the ordinary 
("baryonic") matter, the effective energy density p e ff contains an additional con- 
tribution (87r G ) 1 V 2 A/"(x) = (87r G ) 1 V 2 G(x) due to the position dependence 
of Newton's constant. Since it acts as a source for <f) on exactly the same footing 
as p it mimics the presence of "dark matter" . 

Up to this point the discussion applies to an arbitrary prescribed position 
dependence of Newton's constant, not necessarily related to a RG trajectory. In 
the case of spherical symmetry the natural choice of the geometric cutoff is k = £/r 
with £ a constant of order unity. Hence we obtain the position dependent Newton 
constant G(x) = G(r) as G(r) = G(k = £/r). Writing again G = G [1 +7V], G(k) 
should be such that N <C 1 . 

Let us make a simple model of a spherically symmetric "galaxy" . For an 
arbitrary density profile p = p(r) the solution of the modified Poisson equation 
reads 

r 

0( r )=|dr' ^^l + lM{r) (3.2) 



where M.{r) = 4ir J^dr' r' p(r') is the mass of the ordinary matter contained in 
a ball of radius r. On circular orbits test particles in the potential 1)3. 2JI have the 
velocity v 2 (r) = r <p'(r) so that we obtain the rotation curve 

9 . , GM(r) 1 d , .. , 
v 2 {r) = — ±l + - r -U{r). (3.3) 

We identify p with the density of the ordinary luminous matter and model the 
luminous core of the galaxy by a ball of radius ro . The mass of the ordinary matter 
contained in the core is M(r ) = A4q, the "bare" total mass of the galaxy. Since, 
by assumption, p = and hence M.{r) = M.$ for r > ro, the potential outside the 
core, in the halo, is <f>(r) = -G Mo/r +jV(r)/2. 



6 



Mass Profiles and Shapes of Cosmological Structures 



As an example, let us adopt the power law G(k) oc k~ q with q > which 
was motivated in (Reuter 2004a, b). We assume that this fc-dependence starts 
inside the core of the galaxy so that G(r) oc r q everywhere in the halo. For the 
modified Newtonian limit to be realized, the position dependence of G must be 
weak. Therefore we shall tentatively assume that the exponent q is very small. 
Expanding to first order in q we obtain Af(r) = q ln(r). In the halo, this leads 
to a logarithmic modification of Newton's potential: <fi(r) = —GA4o/r + | ln(r). 
The corresponding rotation curve is v 2 (r) = G Ai^/r + q/2. At large distances the 
velocity approaches a constant i>oo = \fqj2 . Obviously the rotation curve implied 
by the fc _9 -trajectory does indeed become flat at large distances - very much like 
those we observe in Nature. Typical measured values of Woo range from 100 to 
300 km/sec, implying q « 10 -6 which is indeed very small. Including the core 
region, the complete rotation curve reads v 2 (r) — G M. (r)/r + q/2. For a realistic 
M(r) its r-dependence is in rough qualitative agreement with the observations. 

Our v 2 (r) is identical to the one obtained from standard Newtonian gravity by 
postulating dark matter with a density pdm oc 1/r 2 . We see that if G(k) oc k~ q 
with q ss 10 -6 no dark matter is needed. The resulting position dependence of 
G leads to an effective density p c g — p + qj (8n Gr 2 ) where the l/r 2 -term, well 
known to be the source of a logarithmic potential, is present as an automatic 
consequence of the RG improved gravitational dynamics. 

Thus it seems that if the observed non-Keplerian rotation curves are due to 
a renormalization effect, the scale dependence of Newton's constant should be 
roughly similar to G(k) oc k~ q . Knowing this, it will be the main challenge for 
future work to see whether a corresponding RG trajectory is actually predicted by 
the flow equations of QEG. For the time being an ab initio calculation of this kind, 
while well-defined conceptually, is still considerably beyond the state of the art as 
far as the technology of practical RG calculations is concerned. In performing such 
calculations it might help to rewrite the nonlocal terms generated during the flow in 
terms of local field monomials by introducing extra fields besides the metric. This 
is a standard procedure in the Wilsonian approach which often allows for a simple 
local description of the effective IR dynamics. It is tempting to speculate that 
the resulting local effective field theory might be related to the generalized gravity 
theory in (Papers I) which includes a Kalb-Ramond field; it is fully relativistic 
and explains the galaxy and cluster data with remarkable precision. 
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